m1 <- lm(violence ~ days_scaled * treatment, stan_violence, weights = weight)
summary(m1)
confint(m1)

m1a <- lm(violence ~ days_scaled * treatment+age+race6+educ4+gender, stan_violence, weights = weight)
summary(m1a)

m2 <- lm(norms ~ days_scaled * treatment, stan_norms, weights = weight)
summary(m2)
confint(m2)

m2a <- lm(norms ~ days_scaled * treatment+age+race6+educ4+gender, stan_norms, weights = weight)
summary(m2a)

m3 <- lm(affpol ~ days_scaled * treatment, stan_affpol, weights = weight)
summary(m3)
confint(m3)

m3a <- lm(affpol ~ days_scaled * treatment+age+race6+educ4+gender, stan_affpol, weights = weight)
summary(m3a)

#equivalence_test(m1)
#equivalence_test(m2)
#equivalence_test(m3)

stargazer::stargazer(m1,m1a,m2,m2a,m3,m3a, type = 'latex',
                     title="Main Effects",
                     omit.stat = c("ser", "adj.rsq", 'rsq'), 
                     font.size = "scriptsize", column.sep.width = "6pt",
                     df = F, no.space = T,
                     dep.var.labels = c(
                       "Violence", 
                       "Norms", 
                       "Aff Pol"))